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Abstract 

The time evolution of the two-time conditional probability of the classical 
stochastic process is described in an analogous form of the quantum mechan- 
ical wave equations. By using it, we emulate the same strange behaviors as 
those of the weak value in the quantum mechanics. A negative probability 
and abnormal expectations of some quantities remarkablely larger than their 
inherent norms are found in an example of a stochastic Ising spin system. 

Keywords: Weak value, Stochastic process. Two-time conditional 
probability. Stochastic Ising model 

1. Introduction 

A notion of the weak value proposed by Aharonov et al [H has brought a 
new understanding on the quantum observation, i.e. a weak measurement jif 
which hardly disturbs the quantum state. The reason of this strange nature of 
the quantum measurement is that the weak value is defined as an expectation 
with the condition of two-time observations of the initial and the final states 
which differ from one another. This condition is very rare case with little 
probability and is far from the main behavior of a given quantum system. 
Then the observation of the weak value does not disturb the quantum system 
not so fatally. As a result of this rather fictitious probability, the weak value 
happens to be abnormally enhanced from its inherent norm. 

The purpose of this letter is to make the mechanism of this abnormal 
behavior clearer by using a classical stochastic model, in which we can avoid 
the ambiguity of the complex probability in the quantum case 
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We introduce a conventional transformation of the stochastic master 
equation to a self- adjoint form in the following section. A good analogy 
with the quantum mechanics is found by applying this tarnsformation to the 
two-time conditional probability (TTCP). This is shown in Sec. 3. An exam- 
ple of the stochastic Ising model which shows an abnormal enhancement of 
the expectations of some quantities with respect to TTCP is given in Sec. 4. 
In Sec. 5 we discuss an extention of TTCP to a density matrix form to com- 
plete the analogy with the quantum mechanics. The last section is devoted 
to brief summary and discussions. 



2. Self-adjoint form of stochastic master equation 



First let us survey the well-known transformation to a self-adjoint form 
of the stochastic master equation. 

Let a; be a set of stochastic variables described by a time-dependent con- 
ditional probability, P{x,t\xi,ti) for t > ti, which obeys the following sta- 
tionary, Markoffian master equation, i.e. the Chapman- Kolmogorov forward 
equation. 
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The matrix L has an eigenvalue Aq = correspondding to the steady state, 

Pq[x) = lim P{x,t\xi,t{). 

Let us introduce a wave function related to this forward conditional proba- 
bility by 

^(a;, t\x,, U) = <Pq{x)-^P{x, ti), it > U) (2) 
where (I)q{x) = Pq{xY^'^. This function obeys the forward wave equation. 
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where H is defined by 

H(x,x') 



x) L{x, x')(f)o{x') 
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For the time being the initial condition {xi,ti) in is abbreviated. The 
function (f)o{x) is an eigenfunction of Eq.([3]) for Aq = 0. 

The merit of this transformation is that the eigenvalue problem of a given 
master equation is simplified, if the matrix H is symmetric, i.e. 

H{x,x') = H{x',x). 

This situation is widely expected when the detailed balance codition, i.e. the 
time-reversal symmetry |5|, 

Po{x)W{x x') = Po{x')W{x' x), 



or equivalently. 



L{x,x')Po{x') = L{x',x)Poix) 



(5) 



is satisfied. In this case the eigenvalues of H are all real, and non-negative, 
if the steady state is stable. Therefore, (po{x) is the ground state. 

A useful example is the Fokker-Planck equation for a single, continuous 
stochastic variable x. 
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which describes a one-dimensional Brownian motion in a potential F{x) with 
a small diffusion constant e. By using its steady state solutions, 

Po{x) oc exp [—2F{x)/e] and 4>o{x) oc exp [— F(x)/e] , 

we find the continuous variable version of the above formulations, 
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Thus the Fokker-Planck equation is transformed into a self-adjoint form of 
an imaginary-time Schrodinger equation. 



' 2 dx' 
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Figure 1: Stochastic decay process of the metastable state. 



and its eigenvalue problem results in a familiar one of the quantum mechanics. 

Figure. 1 shows an early application \d\ to the so-called Kramers escape 
problem [7|. The stochastic decay rate of the metastable state in a double- 
well potential F{x) is given by the first excited eigenvalue Ai of the cor- 
responding Schrodinger potential V{x). The first excited state is almost 
degenerated with the ground state for the small diffusion constant e. 

3. Two-time conditional probability 

So far the quantum mechanical reformulation merely helps us to simplify 
the eigenvalue problem of a given master equation. None of the remarkable 
quantum mechanical phenomena appears, until we are concerned with the 
TTCP, 

P{x,t\xi,ti; X{,tf), ti < t < t{ . ( ; denoting 'and') (8) 

By using the Markoffian property and the well-known equality of the simplest 
Bayes' theorem, 

P{A nB) = P{A\B)P{B) = P{B\A)P{A), 

repeatedly, the TTCP can be written in the following form with a pair of the 
wave functions as 

P{x, t\xi, U; X{, tf) = ] ^(a;, t\x{, tf)4j{x, t\xi, U), (9) 

where the conjugate wave function ijj is related to the so-called posterior con- 
ditional probability, P{x,t\xi,tf) for t < tf, by 

i>{x, t\xi, tf) = (j)o{xy^p{x,t\xi, tf), (10) 
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and obeys the backward wave equation, 

^Hx,t) = J2H\x,x')^ix',t). (11) 
x' 

Here if^ is the hermite conjugate of H, i.e. the transposed matrix in the 
present case. The eigensystem is common with the forward equation Eq.(l3]), 
when H is hermite, i.e. real and symmetric as has been assumed here. 

The denominator in Eq. (Q is the weight of overlap between the two wave 
functions defined by an inner product, 

(^fl^i) = ^^(a;',t|ajf,tf)7/^(a;',t|a;i,ti)- (12) 
x' 

Of course this quantity is real, while the corresponding quantity in the quan- 
tum mechanics is complex in general. 

Let us define the ket- and the bra-vectors by 

iV'i(^)) = {i^ix,t\xi,ti)}^ and (^f(t)| = {^p{x,t\xi,ti)}. 

Then the wave equations Eqs.© and f|TT]) are rewritten in the quantum 
mechanical form as 

l\^,{t)) = -H\Ut)) and l{^,{t)\ = {Mm. (13) 

Henceforth, H is called the Hamitonian. 

By using this pair of the Schrodinger equations it is shown that the overlap 
integral, {i'di'i) given by Eq. (fT^ does not depend on the current time t, i.e. 

= (Mmmt)) - mmim) = o. 

Further it can be shown that this overlap integral has the following properties 
in the two limits; 

(i) lim (^fl^i) = 1, 

(14) 

(ii) Mm (ipfliJi) = [(j)o{xf)(l)o{xi)] ^6{xi-Xi). 
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Note that the two-time conditional expectation (TTCE) of a physical 
quantity Q with respect to TTCP defined by 



(i;f) 



^ Q{x')P{x\ t\x;, U; Xi, ti) 



X 



(15) 



has just the analogous form of the weak value in the quantum mechanics jsj. 

Thus the TTCP is a non-linear quantity composed of a product of a pair 
of the forward and the backward wave functions, and cannot be described by 
a closed, linear evolution equation. Then it happens that the principle of the 
probability superposition is violated and the interference of wave functions 
may occur. However, its example is omitted here because none of nontrivial 
phenomenon from this view point has been found, yet. The reason may be 
that the wave functions are always real and possitive in the present case. Let 
us discuss only the weak value in the rest. 



4. Stochastic model of classical Ising spins 

An example is a pair of the classical Ising spin cr = ± 1 having an exchange 
interaction, 

E{x) = -J(TiO-2, 

where x = ((Ti,cr2)- Let us number the stochastic variable x in the order, 
(1, 1), (1, —1), (—1, 1), (—1, —1) and choose the following transition matrices, 
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where p = e' 



Evidently this transition matrix W satisfies the detailed 



balance condition. 



at the steady state, i.e. the thermal equilibrium of a temperature parameter, 
/3 = 1/ k-QT. With use of the equilibrium distribution function. 



we find the corresponding hermite Hamiltonian, 
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= (1 + p^) (7o (g) (Jo - (1 - p^) (g) (72 - p ((7o (g) (7j; + (7a; (8) CTq), (17) 

where (Xa. and are the usual Pauli matrices and (7o denotes the two dimen- 
sional unit matrix 72- This is the Hamiltonian of a pair of quantum Ising 
spins with the exchange interaction in a transverse magnetic field. 
The eigenvalues and the eigenstates of this Hamiltonian 
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can be easily obtained, where |0) = |</)o), the ground state. Here the famihar 
notations t, i ai'c used for a = ±1. Note that the first excited state is almost 
degenerated with the ground state for a small transition probability p"^. 

By using this eigensystem we can calculate the state vectors, iV'i(^)) and 
{ij^i{t)\ for arbitrary initial and final states in just the same manner of the 
elementary quantum mechanics except for the fact that the time t is imagi- 
nary. 

An interesting example is the case where the initial and the final states 
differ from each other, just like the case of the weak value. For example, let 



that is. 



= (tt) at t = and Xf ^ {H) at t ^ tf, 
P(a;,0) = (1,0,0,0) and P(a;, if) = (0, 0, 0, 1), 
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Figure 2: Two-time conditional probability 



or equivalent ly, 

|^i(0)) = v/2(l+p2) Itt) and {Mt()\ = + p') (Ul- 

By using the eigenvector expansion we obtain, 

\Ht)) = |0) + VTTp2 e-^i* |1) +Pe-^'* |3), 



;i8) 



(^f(t)l = (0| - v^TTp^ e-^^M(l| +p e-^3(*f-*)(3|, 

and 

(V'flV'i) = 1 - (1 + p^) e-^i*f + p2 e-^^*f (> 0). (19) 

The TTCP is shown in Figure. 2. This result itself is very natural and 
well-expected, all probabilities being always non-negative. 

A strange behavior appears when we use the basis {\k), k = 0,1, 2, 3}, the 
eigenstates of the Hamiltonian H instead of the spin states {\x) = |cti(T2)}. 
We can calculate the probability, i.e. the TTCE of the projection operator 
onto each eigenstate k in the same manner. The result is given by 



P(0,t) 



(^K^)|o)(o|^i(t)) ^ 1 



P{2,t) 
P(3,t) 



(^fl^i) (^flV'i) 

(^fi^i) 

(^f(t)|3)(3|^i(t)) ^ p^e-^3tf 



(20) 
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Figure 3: Abnormal and normal TTCE of for p = 0.2 and p^tf = 0.01. 



The fictitious negative probability is found in P(l,t). Of course the com- 
pleteness of the probability, 



3 

E 

fc=0 



is satisfied evidently because of Eq. (fT9|) . 

A related unusual behavior to this fact is the abnormal enhancement of 
some observables. For example, if we calculate the TTCE of a quantity, 

= ^(o-x ® o-Q + ao ® 0-^), (21) 

an abnormal behavior 



(V^flV^i) 



> 1, 



is found for sufficiently small p and t{. An example is shown in Figure. 3. Note 
that the natural norm of must be less than 1, because the eigenvalues of 
Mj, are { — 1,0,0,1}. When the transition rate is very small, i.e. p'^tf <^ 1, 
we find 

(M,)(i.f) > 1. 

A plain reason of this singular behavior is that the overlap integral in 
the denominator may be expected to be very small owing to (ii) of Eq.f ll4p . 
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Figure 4: Abnormal and normal TTCE A = ax <^ (Tx for p = 0.2. 



whenever the initial and the final states differ from each other, i.e. x-^ ^ Xf. 
This means that to reach Xf = starting from X[ = in a given time 
occurs scarcely and is far from the main flow of the conditional probability. 
On the contrary none of such strange behaviors are found when Xi = Xf, e.g. 
Xi = Xf = ('I't) • The result for the latter case for the same parameters as the 
upper abnormal case is shown by the lower curve in Figure. 3, its maximum 
being ~ 0.09 at t = tf/2 and minimum ~ 0.05 at t = and t{. 

In Figure. 4 the TTCE of another quantity A = ax ^ are shown also. 
Note that A is commutative with H and a conserved quantity. Then the 
horizontal axis in this figure shows a parameter of the transition probability, 
not the time. 

5. Extension of TTCP to a density matrix 

It should be noted that the physical quantities and A are non- diagonal 
in the spin-state representation and have no corresponding quantities in the 
classical Ising spin system. They are related to the transition rate of the 
stochastic Ising spin. In order to calculate the expectations of such non- 
diagonal quantities we need an extension of the TTCP to the two-time con- 
ditional (TTC) density matrix defined by 



P(i;f)(^) 



1 



\m){Mt)\ 



(^fi^i) 
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x,x' 

From the definition Eq. (fT2|) of tlie overlap integral {ipi\il)-^, it is evident that 

Tr P(i;f)(t) = / , ! , \ V^(a;, t|a;f, tf)^(x, 0) = 1. 
(V'fl^i) ^ 

It should be noted, however, that diagonal elements of this density matrix are 
not always positive as is shown by Eq. fl^U]) in Sec. 4, when it is diagonalized 
by using the basis = 0, 1, 2, 3}, the eigenstates of the Hamiltonian H. 

With use of this density matrix the definition Eq. flT^ of the TTCE is 
extended as 

(Q)(i;f) = Tr P(i;f)(5- 

This definition of the TTCE results in the classical one, if Q is diagonal. 

The notion of this density matrix has not been used in the conventional 
classical stochastic process. It should be emphasized, however, that this 
quantity is within a scope of the classical stochastic process itself, because 
the wave functions, ip and in Eq. fl22l) are related to the forward and the 
posterior, classical conditional probabilities, respectively. In addition, we 
have an alternative expression for ■?/', 

lp{x' ,t\xuU) = ij{x' M^ut) {=M^')-^P{x',U-t\xu^)), (23) 

or equivalently, 

P{x' ,t\xuU)PQ{x,) = P{xuU\x' ,t)Po{x') 

= P{x',U\xut)P^{x,), (24) 

for t < t{ due to the time-reversal symmetry corresponding to the detailed 
balance. Then the density matrix Eq. fl22|) can be written as 



1 ^ P{x',t{-t\x{,0)P{x,t\xi,0) 
P(i;f) W = , , I , , > , . M , . X \x){x'\. 25 



x,x 



This fact means that we can define the TTC density matrix with only a 
pair of the usual, forward conditional probabilities for two individual initial 
states, Xi and Xf. 
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6. Summary and discussions 

Except for the facts that the time is imaginary and the wave function 
is always real and positive, the classical stochastic process can be described 
in an analogous form of the quantum mechanics, if we use the TTCP. For 
example, the abnormal behaviors of the weak value in the quantum mechanics 
are emulated. Note that the TTCP and the weak value are always real in the 
present classical case. Therefore, the origin of such abnormal behaviors is 
clearer than the quantum mechanical case where complex quantities appear. 
In addition, if we have not the explicit solution of the eigenvalue problem, we 
may calculate the weak value at least with use of the Monte-Carlo simulation 
which is often used to investigate the stochastic model. 

The importance of the weak value in the quantum mechanics is that it is 
related to the new notion of the weak measurement without disturbing the 
quantumu state. An analogous notion of the latter in the classical stochastic 
process has not been found yet. 
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